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Effects of spontaneous toroidal ordering on magnetic excitation are theoretically investigated for a localized spin 
model that includes a staggered Dzyaloshinsky-Moriya interaction and anisotropic exchange interactions, which arise 
from the antisymmetric spin-orbit coupling and the multi-orbital correlation effect. We show that the model exhibits 
a Neel-type antiferromagnetic order, which simultaneously accompanies a ferroic toroidal order. We find that the oc¬ 
currence of toroidal order modulates the magnon dispersion in an asymmetric way with respect to the wave number: a 
toroidal dipole order on the zigzag chain leads to a band-bottom shift, while a toroidal octupole order on the honeycomb 
lattice gives rise to a valley splitting. These asymmetric magnon excitations could be a source of unusual magnetic re¬ 
sponses, such as nonreciprocal magnon transport. A variety of modulations are discussed while changing the lattice and 
magnetic symmetries. Implications of candidate materials for asymmetric magnon excitations are presented. 


Effects of the spin-orbit coupling (SOC) on magnetic mate¬ 
rials have attracted a great interest for the rich physics emer¬ 
gent from entangled spin and orbital degrees of freedom. 
Magnetoelectric effects caused by coexistence of ferromag¬ 
netic (FM) and ferroelectric orders are typical examples. 
Another examples are magnetotransport properties, such as 
the topological Hall effect in magnetic skyrmions."^ In these 
phenomena, the lack of spatial inversion symmetry plays an 
important role: the SOC acquires an antisymmetric compo¬ 
nent with respect to the wave number, which is called the an¬ 
tisymmetric spin-orbit coupling (ASOC). 

The ASOC can exist even in centrosymmetric systems in a 
hidden form. Suppose the inversion center is located at an off¬ 
site position, and the inversion symmetry is broken at each lat¬ 
tice site. Then, the local asymmetry may induce an ASOC in 
a sublattice-dependent form, whose net component vanishes.^ 
For instance, this class of centrosymmetric systems are found 
in zigzag chain, honeycomb, and diamond lattices. 

In such situations, interesting physics arises once a spon¬ 
taneous breaking of sublattice symmetry takes place by elec¬ 
tron correlations: a net component of ASOC is induced by the 
additional inversion symmetry breaking.This opens the 
possibility of controlling the ASOC through the electronic de¬ 
grees of freedom. Furthermore, such an emergent ASOC has 
a close relationship to odd-parity multipoles. For instance, 
in the one-dimensional (ID) zigzag-chain system where the 
ASOC is present in the staggered form, a Neel-type antiferro¬ 
magnetic (AFM) order induces an additional net ASOC. This 
accompanies a ferroic order of odd-parity multipoles, such as 
magnetic toroidal and quadrupole, which break both spatial 
and time-reversal symmetries.^’^ Interestingly, such a toroidal 
order modifies the electronic band structure in an asymmet¬ 
ric way in the momentum space, resulting in unusual off- 
diagonal responses, such as the magneto-current effect.^’^ 

Spontaneous toroidal ordering in the presence of ASOC is 
also expected to affect dynamical properties, such as a non¬ 
reciprocal directional dichroism, through the collective exci- 
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Eig. 1. (Color online) Schematic pictures of the zigzag chain along the x 
direction for (a) the localized spin model in Eq. (1) and (b) the itinerant elec¬ 
tron model in Eq. (2). (c) Schematic picture of the honeycomb lattice for the 
localized spin model in the xy plane. In (a) and (c), A and B represent two 
sublattices. 


tations. Moreover, as the odd-parity multipoles are associated 
with conventional orders, such as the Neel AFM, they are re¬ 
alized without magnetic frustration that usually suppresses the 
relevant energy scale. Hence, they have a potential for next- 
generation electronic devices working at higher temperatures. 
Nevertheless, the microscopic theory has not been fully elab¬ 
orated thus far. 

In the present paper, we investigate the magnon excita¬ 
tions in Neel-type AFM ordered states in the centrosymmetric 
systems with local asymmetry. Considering the Heisenberg 
model with a staggered Dzyaloshinsky-Moriya (DM) inter- 
action^^"^"^ and anisotropic exchange interactions,^^’which 
originate from the ASOC and multi-orbital correlation effect, 
we show that the AFM orders with ferroic toroidal ordering 
can be realized. We find that the spontaneous odd-parity mul¬ 
tipole orders result in peculiar asymmetric deformations of 
the magnon dispersion, such as a shift of the band bottom and 
a valley splitting. Such deformations lead to unconventional 
off-diagonal responses even in the absence of external fields. 
Extending the analysis to other cases including FM and lo¬ 
cally symmetric cases, we systematically elaborate how the 
magnon excitation is modulated by the lattice and magnetic 
symmetries. 

Let us begin with the ID zigzag chain, the simplest realiza- 
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tion of the centrosymmetric lattices with local asymmtry. We 
consider the Heisenberg model supplemented by anisotropic 
exchange interactions and a sublattice-dependent DM interac¬ 
tion, whose Hamiltonian is given by 

« = 2 J„S, ■ Si + 2G,J(S-S; -sfs; - s’s;) 

ij ij 

+ \(Sap X Sap+i) - {Sbp XS'bp+i)] • z, (1) 
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where Si is the 5 = 1/2 operator at site i = (/, p) (I and p 
denote the indices for the sublattice A and B, and unit cell, re¬ 
spectively). The first term represents the Heisenberg-type ex¬ 
change coupling for nearest-neighbor (NN) and next-nearest- 
neighbor (NNN) spins, Ji and J 2 , respectively. The second 
term in Eq. (1) is the anisotropic exchange interactions^^’ 
for NN and NNN spins, Gi and G 2 , respectively. The former 
Gi originates from the multi-orbital correlation effect, while 
the latter G 2 from the ASOC, as discussed below. The third 
term in Eq. (1) is the antisymmetric exchange interaction be¬ 
tween NNN spins (DM interaction)[see Eig. 1(a)], whose 
origin is also the ASOC as shown below. Note that there is no 
DM interaction for NN spins because of the inversion center 
at the bond median. 

In the second and third terms in Eq. (1), the sublattice- 
dependent ASOC inherited from the local asymmetry of the 
zigzag lattice plays an important role. This is understood, for 
instance, by considering the strong-coupling limit of an effec¬ 
tive single-band model,^ whose Hamiltonian is given by 

= E + H.c.) + uYj 

ij,cr i 

+ E Sik) • 0-cr<r' • (2) 

k,(T,(T' 

The first and second terms comprise the Hubbard model in a 
standard notation; ti and t 2 are NN and NNN hoppings, re¬ 
spectively. The third term represents the sublattice-dependent 
ASOC in the form of the spin-dependent electron hopping, 
which originates from the atomic SOC, local asymmetry, and 
multi-orbital hybridization effects between orbitals with dif¬ 
ferent parity .^’g is the asymmetric vector with respect 
to k: g(k) = 2asin(ka)z [a is the amplitude, a is the lat¬ 
tice constant (a = 1), and z is the unit vector along the z 
axis; see Eigs. 1(a) and 1(b)], and a is the Pauli matrix. The 
sign factor comes from the zigzag structure, as an electron at 
the A and B sublattice sites feels an opposite potential gradi¬ 
ent. By considering the strong-coupling limit of the model in 
Eq. (2) at half-filling, we obtain the coupling constants in the 
model in Eq. (1) by the second-order perturbation in terms 
of tju (m = 1,2) and ajU as G 2 = 4a^lU, 

and D = 2 On the other hand, in general, the multi- 

orbital correlation effect, which is not taken into account in 
Eq. (2), favors the Neel-type AEM ordering along the z direc¬ 
tion in the presence of the atomic SOC and the multi-orbital 
hybridization, as demonstrated in sl d-p model^^ and an ex¬ 
tended Kondo lattice model.^ Such a multi-orbital correlation 
effect is incorporated in the Gi term in Eq. (1). Hereafter, we 
consider /i > 0 and J 2 > 0 with taking Ji as an energy unit, 
and Gi > 0, G 2 > 0, and D = 2 V/ 2 G 2 . 

Eirst, we examine the ground state of the model in Eq. (1). 
To stabilize the long-range orders, we here consider a three 


Eig. 2. (Color online) The ground-state phase diagrams for the model in 
Eq. (1) at (a) /2 = 0 and (b) J 2 = 0.2. The contour displays the ordering wave 
number q *. xy(z)-AFM represents the AFM state with the spin polarization 
along the xy{z) direction. 


dimensional system composed of ferromagnetically weakly- 
coupled ID zigzag chain. We obtain the ground-state phase 
diagram by the Luttinger-Tisza method by treating the spins 
as classical vectors.The results are plotted in Eigs. 2(a) and 
2(b) for /2 = 0 and J 2 = 0.2, respectively. The contour dis¬ 
plays the optimal wave number (f. When J 2 = Gi = G 2 = 0 
[the origin in Eig. 2(a)], the stable spin configuration is a ^ = 0 
collinear AEM without any spin anisotropy. With increasing 
Gi (G 2 ), the NN (NNN) anisotropic interaction stabilizes the 
AEM ordering along the z direction (z-AEM) [in the vy plane 
(vy-AEM)] without any spin canting. Note that the system has 
the spin-rotational symmetry around the z axis. With further 
increasing G\ and G 2 , incommensurate magnetic orders with 
longer period ^ A 0 appear. As shown in Eig. 2(b), the ^ A 0 
region becomes larger for larger J 2 , as there is frustration be¬ 
tween Ji and J 2 . 

The results in Eigs. 2(a) and 2(b) show that the model ex¬ 
hibits the collinear z-AEM order for sufficiently large G 1 /G 2 . 
As described above, the z-AEM order on the zigzag chain 
is of particular interest as it accompanies the odd-parity 
multipoles, such as magnetic toroidal dipoles and magnetic 
quadrupoles, through the spin-orbit coupling effect when D 9 ^ 
0. The magnetic toroidal dipoles we are considering in this 
study are defined by ^ where ri is the position 

vector from the inversion center to the lattice site In such 
a situation, the z-AEM order possess the toroidal dipoles in 
the V direction because Si || z and r, || y by taking the bond 
median as the origin. 

Next, we examine magnetic excitations in the collinear 
z-AEM state by using the spin-wave theory. We adopt the 
standard Holstein-Primakoff transformation, which is repre- 
sentedby5+^ = yflSuqiS'^^ = V25Z?J), 5^^ = 

= yl^bq), andS\^ = S - aluq (5^^ = -5 -I- blbq) for 
the A (B) sublattice. Here, a and h are the boson operators for 
the A and B sublattices, respectively, and q is the wave num¬ 
ber along the chain direction. We adopt the linear spin-wave 
approximation, in which the magnon-magnon interactions are 
ignored. Then, the spin-wave Hamiltonian in the q space is 
obtained as 

'H = Eo + 

q 

+ E +4^4)’ (2) 
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Eig. 3. (Color online) (a) Magnon dispersion in Eq. (4) in the z-AFM state 
on the zigzag chain. The result is calculated at (/ 2 , , G 2 ) = (0.1,0.05,0.02). 
For comparison, the results for (/ 2 ,Gi,G 2 ) = (0,0,0) and (/ 2 ,Gi,G 2 ) = 
(0,0.05,0) are shown by the dotted and dashed curves, respectively. All the 
dispersions are doubly degenerate. J 2 -G 2 plots of (b) and (c) As at Gi = 
0.05. The blank regions represent the areas where the xy-AFM and other 
incommensurate orders appear. 


where = 25 [(/i + Gi) +/ 2 (cos^ - 1) - G 2 (cos^ + 1)], 
= ISD^mq, Bq = 2S{Ji - Gi)cos(^/2), and Eq = 
-NS^{J\ - /2 + Gi - G 2 ). By using the Bogoliubov transfor¬ 
mations aq = aq cosh 6q + b^_q sinh 6q and= aq sinh 6q + 
blqCO^hOq [Oq = j 2) qj we can diagonalize 

Eq. (3) into the form of E( = G^qia^a^ + Pl^q) + const., 
where the magnon dispersion relation ojq is obtained as 

+ (4) 

Figure 3(a) shows the magnon dispersion in the z-AFM 
state at /2 = 0.1, Gi = 0.05, and G 2 = 0.02. The magnon ex¬ 
citation spectrum has a peculiar form: in addition to the gap 
opening, the dispersion undergoes an asymmetric deforma¬ 
tion with respect to q, leading to a shift of the band bottom 
from q = 0. This is similar to the electronic band structure 
in the presence of the toroidal ordering.^’^ For comparison, 
we also show the results at /2 = Gi = G 2 = 0 (dotted lines) 
and J 2 = G 2 = 0 (dashed lines). The results indicate that Gi 
opens the excitation gap and D = 2 V/ 2 G 2 brings about the 
asymmetric deformation. 

We analyze the effect of Gi, G 2 , and J 2 more carefully. Fig¬ 
ures 3(b) and 3(c) represent the band-bottom shift Aq and the 
spin gap As on the J 2 -G 2 planes at Gi = 0.05. The results 
show that J 2 and G 2 significantly affect the asymmetric de¬ 
formation of the magnon dispersion, as shown in Fig. 3(b). 
This is because the deformation is caused by A^^ oc D sin ^ 
originating from the sublattice-dependent DM interaction in 
the third term in Eq. (1). On the other hand, the spin gap de¬ 
pends on both Gi and G 2 , while it does not show strong J 2 
dependence, except in the region close to an incommensurate 
phase, as shown in Fig. 3(c). 

For further elucidating the relation between the magnon 
dispersion and the symmetry of the system, we extend our 
analysis to other cases. Namely, we consider both AFM and 


Eig. 4. (Color online) Summary of the characteristic of the magnon disper¬ 
sions on the ID systems. The row represents staggered- or uniform-type DM 
interaction, while the column represents staggered- or uniform-type mag¬ 
netic ordering polarized in the z direction. The dispersions are calculated for 
(/i,/ 2 ,Gi,G 2 ) = (1,0.1,0,05,0.02) (upper left), (1,0.1,0.05,0.02) (lower 
left), (-1,-0.!,-0.05,-0.02) (upper right), and (0,-1,0,-0.02) (lower 
right) in the model in Eq. (1).^^ While the lower left is obtained by assuming 
the same sign for D (the uniform DM) in Eq. (1), the upper right is calculated 
in the EM ordered state. A schematic picture for each magnetic state is shown 
in each inset. The vertical dotted lines represent the shifted band bottoms. 


FM cases while changing the parameters in Eq. (1) beyond 
the range expected from the perturbation theory. We also con¬ 
sider a straight ID chain in an applied electric field, which has 
a uniform DM interaction. 

Figure 4 summarizes the characteristic of the magnon dis¬ 
persions in the matrix form for the staggered/uniform DM in¬ 
teractions (row) and the AFM/FM orders (column). The above 
result for the spontaneous z-AFM ordering on the zigzag 
chain is shown in the upper-left panel. A similar asymmet¬ 
ric dispersion also appears for the z-FM ordering under the 
uniform DM interaction (lower-right).^^ Note that the latter is 
in the same category as the toroidal magnon discussed for the 
helical magnetic structure in an applied magnetic field.^^’^"^ 
These asymmetric magnon dispersions lead to peculiar mag¬ 
netooptical phenomena, such as the nonreciprocal directional 
dichroism.^^’^^ 

On the other hand, when the z-AFM order occurs on the 
straight chain with the uniform DM, the magnon dispersion 
undergoes a different modulation (lower-left), which is simi¬ 
lar to a “Rashba-type” splitting in the electronic band struc¬ 
ture under a uniform ASOC. Meanwhile, in the case of the 
z-FM order on the zigzag chain, there is a symmetric modula¬ 
tion, reflecting the spatial inversion symmetry (upper-right). 

Finally, let us extend the analysis to the two-dimensional 
case. We study the magnon dispersion in the collinear z-AFM 
state on the honeycomb lattice. Similar to the ID zigzag case, 
the spin model includes the anisotropic and antisymmetric in¬ 
teractions, as shown in Fig. 1(c). D is positive (negative) in 
the counterclockwise (clockwise) direction in the hexagonal 
plaquettes. In the case of the honeycomb lattice, a collinear z- 
AFM order accompanies ferroic toroidal octupoles instead of 
dipoles, reflecting the presence of threefold rotational sym- 
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Eig. 5. (Color online) (a) The ground-state phase diagram on the hon¬ 
eycomb lattice at /2 = 0.1. The contour displays the ordering wave 

number q* = + qy . (b) Magnon dispersion under the z-AFM order¬ 

ing on the honeycomb lattice. The result is calculated at (/ 2 ,Gi,G 2 ) = 
(0.1,0.015,0.005). The dotted and dashed curves represent the magnon dis¬ 
persions at (J 2 ,Gi,G 2 ) = (0,0,0) and (J 2 ,Gi,G 2 ) = (0.1,0.015,0), respec¬ 
tively. The inset of (b) shows the Brillouin zone. 


metry of the lattice, i.e., the net toroidal-dipole component 
is zero. We find that, by the Luttinger-Tisza analysis, such a 
z-AFM order with toroidal octupoles becomes stable in the 
ground state for large G 1 /G 2 , although the region becomes 
narrower than that for the zigzag chain, as shown in Fig. 5(a). 

In the collinear z-AFM ordered state, the magnon disper¬ 
sion does not show any shift of the band bottom, as the lowest 
contribution to the ASOC is in the third order in (qx, qy) due 
to the threefold rotational symmetry. Instead, the asymmetry 
appears near the Brillouin zone boundary around the K and 
K’ points, as shown in Fig. 5(b). This is a valley splitting, 
similar to the electronic band structure discussed in noncen- 
trosymmetric compounds, such as monolayer transition metal 
dichalcogenides.^^ Note that the band-bottom shift is obtained 
also on the honeycomb lattice once a uniaxial pressure is ap¬ 
plied, since it breaks the threefold rotational symmetry and 
induces a ^-linear contribution to the ASOC as in the zigzag 
chain case. 

To summarize, we have clarified that spontaneous order¬ 
ing of toroidal multipoles modulates the magnon excitations 
in an asymmetric way in the momentum space. The obser¬ 
vation of asymmetric magnon spectra will provide an exper¬ 
imental probe for toroidal multipoles. There are many ex¬ 
perimental candidates. One is the AFM zigzag compound 
a-Cu 2 V 207 .^^ Since the compound may possess both the 
staggered and uniform DM interactions due to the peculiar 
magnetic and lattice structures, both a band-bottom shift and 
“Rashba-type” splitting are expected to be observed. An¬ 
other candidates are transition metal tricalcogenides, such as 
MnPSs and MnPSes:^^"^^ the z-AFM order on the honeycomb 
lattice will lead to a valley splitting in the magnon spectrum. 
A distorted honeycomb compound yS-YbAlB 4 is also a can¬ 
didate, as it shows the AFM ordering under a uniaxial pres¬ 
sure a band-bottom shift is expected in addition to a valley¬ 
splitting. We note that the diamond-lattice systems will also 
be in the present scope, as they possess the similar physics 
related with the local asymmetry. 

Similar asymmetric excitations will be obtained also in itin¬ 
erant electron systems, such as the Hubbard and Kondo lattice 
models, when the sublattice-dependent ASOC is present. Re¬ 
flecting the itinerant electron degree of freedom, asymmet¬ 


ric Stoner-type excitations are expected in addition to asym¬ 
metric magnon excitations. Such an analysis is left for future 
study. 
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